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ABSTRACT 

We present the results of applying a percolation algorithm to the initial release of the Two Micron 
All-Sky Survey Extended Source Catalog, using subsequently measured redshifts for almost all of the 
galaxies with K < 11.25 mag. This group catalog is based on the first near-IR all-sky flux-limited 
survey that is complete to \b\ = 5°. We explore the dependence of the clustering on the length 
and velocity scales involved. The paper describes a group catalog, complete to a limiting redshift of 
10 4 km s _1 , created by maximizing the number of groups containing 3 or more members. A second 
catalog is also presented, created by requiring a minimum density contrast of 5p/p > 80 to identify 
groups. We identify known nearby clusters in the catalogs and contrast the groups identified in the 
two catalogs. We examine and compare the properties of the determined groups and verify that the 
results are consistent with the UZC-SSRS2 and northern CfA redshift survey group catalogs. The 
all-sky nature of the catalog will allow the development of a flow-field model based on the density 
field inferred from the estimated cluster masses. 

Subject headings: galaxies: clusters: general — galaxies: distances and redshifts — large-scale struc- 
ture of universe 



1. INTRODUCTION 

The Two Micron All-Sky Survey (2MASS) began in 
the early 90s with the purpose of mapping the Milky 
Way and nearby Universe. Previous all-sky surveys suf- 
fered from a variety of selection effects, many of them 
inconsistent across the sky. Observations at optical 
wavelengths suffer from severe extinction at low galac- 
tic latitudes, motivating work using surveys conducted 
in the infrared. Many IRAS-selected galaxy samples 
have been in vestigated as tracers of the galaxy den - 
sity field, e.g. IStraniss et a.U Ip92^ : IPisher etaTI ll 995). 
bas ed on the 1.9 Jy and 1.2 Jy samples, respectively, 
and IBran chini et alJ IU999T ) based on the PSCz catalog 
ijSaunders et alJfeOOOl) . however these samples are based 
on fluxes in the far-infrared and miss many early-type 
galaxies, thus underestimate the total galaxy number 
density. Even though the IRAS-selected samples are not 
biased by extinction, they still suffer from confusion in 
high-density regions. 

Since galaxies' spectra peak at ~1.6 /im, a survey in 
the near-infrared is optimized for their detection at the 
flux limit of the survey. By sampling uniformly over the 
entire celestial sphere in the J (1.25 /xm), H (1.65 /im) 
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and K (2.16 ^m) bands, 2MASS has been designed to 
maximize the number of galaxies detected at a specified 
flux limit, producing the most complete all-sky survey 
performed to date r2MASS. ISkrutskie et aLll2006tl . 

Two catalogs, complete to K ~ 13.5 mag, were re- 
lease d in early 2003 (se e th e 2MASS explanato ry supple- 
ment iCutri et alJl2003t and lJarrett et allfe OOO): a point- 
source catalog with 470,992,970 entries and an extended 
source catalog (XSC) with 1,647,559 objects classified as 
galaxies. Although designed for completeness down to 
low galactic latitudes, the 2MASS XSC still suffers from 
confusion near the galac tic plane. The 2MASS Redshift 
Survey (hereafter 2MRS, Huchr a et alJl2005blaf) uses the 
XSC as its input master list and aims to produce an 
all-sky, (extinction-corrected) flux-limited redshift cata- 
log that will eventually be complete to K = 13.0 mag 
above \b\ = 5°. The 2MRS is currently 99.9% complete 
to K = 11.25 mag and \b\ > 5°. 2 

In this paper, we create a redshift-limited catalog of 
groups, uniformly sampled from the entire sky. By as- 
suming the identified groups are virialized systems, we 
are able to provide estimates of the group masses, avoid- 
ing the necessity to assume an intrinsic mass-to-light ra- 
tio. The local Universe is sufficiently inhomogeneous at 
the scales in question that the dynamics due to our in- 
teractions with nearby groups are non-negligible. Due to 

2 See 2MRS data release, Huchra et al. (2006), In preparation. 
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the nature of the all-sky group catalog presented here, we 
will now be able to estimate the local density field due 
to baryonic matter in the local Universe. A flow-field 
model produced from this catalog can be used in con- 
junction with observations in order to answer the ques- 
tion of whether baryonic matter is a genuine tracer of 
dark matter. 

The creation of group catalogs is not a new concept, 
however the methods employed in developing these cat- 
alogs have evolved with the enhancements in instru- 
mentation. Early group ca talogs were based on lim- 
ited or subjective data (e.g Ide Vaucouleura 11975)) . as- 
sociating members based on similarity in apparent mag- 
nitude, positional coin cidence and (if available) redshift. 
iTurner fc Gottl 1)197 6') proposed a method that identifies 
regions in which the surface number density on the sky is 
enhanced, creating group catalogs from two-dimensional 
data. This technique suffers because the typical angu- 
lar separation of galaxies in a group will vary with dis- 
tance, thus nearby groups with large angular radii will 
not be identified. Furthermore, when applied to flux- 
limited surveys, this method will identify different groups 
for different limiting fluxes. 

More recently, the use of objective algorithms to 
identify groups based on both their position on the 
sky a nd in redshift space has become widely accepted 
(e.g. iHuchra fc Gelled Il982l iGeller fc Huchral 11983 
IRamella et al.||1997|: i Diaferio et al.lll999t iGiuricin et al.l 
120001 iRamella et all 120021) . using methods designed to 
find the same groups regardless of the limiting magni- 
tude of the sample. The applicability of a particular 
group-finding algorithm depends on the properties of the 
sample in question. For example. iMarinoni et alJ ( 2002) 
show that the Vornoi-Delauney method successfully re- 
produces the distribution of groups in velocity dispersion 
in a mock sample based on the Deep Extrag alactic Evo- 
lution ary Probe (DEEP2) Redshift Survey (iDavis et all 
12003ft : this method is adapted by iGerke et al.l (|2005f) 
for application to the DEEP2 sample. The SPSS team 
developed an algorithm (C4, iMiller et all [2005) that 
searches for groups in th ree space-dimens i ons as well as 
four photometric colors. Kochanck et al.l l)2003|) used a 
matched filter algorithm t o study clus t ers i n 2MRS at the 
89% completeness level. lYang et all 1)20051) have devel- 
oped a halo-based group-fi nder and successfully a pplied 
it to the 2dFGR S sample l)Merchan fe Zandivarezll2002l 
lEke et all 120041 ) The same technique h as been applied 
to the SDSS hv lWeinmann et a.U pfflfiT) . 

In this paper we apply a variable linking-length perco- 
lation (als o commonly referred to as a friends- of- friends) 
algorithm l)Huchra fc Gellerll982l hereafter HG82) to de- 
termine the groups present in 2MRS. The velocity disper- 
sion within the groups will allow estimates of the virial 
masses of the groups, thus providing a method to trace 
the density field associated with luminous matter in the 
local Universe. 

We begin with an outline of the group-identification al- 
gorithm in below. We discuss the modifications made 
to the data sample prior to the application of the al- 
gorithm in SI presents a discussion on the choice 
of parameters used in the group-identification algorithm. 
The group catalogs and their properties are discussed in 
§5\ and we summarize our conclusions and discuss the 
potential applications of the catalogs in §§\ 



2. GROUP-IDENTIFICATION ALGORITHM 

We use the algorithm described in HG82 to identify 
groups of galaxies in the K < 11.25 mag version of the 
2MRS catalog. The procedure is outlined briefly below. 

We compare each galaxy in the catalog with its neigh- 
boring galaxies; for each pair of galaxies a linking length 
-Di(V r av g ) is computed that depends on the average red- 
shift of the galaxies, V avg . Given two galaxies with 
an angular separation, 9, we ask whether their pro- 
jected separation, D\i — sin(9/2)V avg / H , is less than 
-D/^Uavg). If this is true, and the difference in redshift, 
Vu = \Vi — V%\, is less than some linking velocity, Vl, 
then we identify both galaxies with the same group. Dl 
is defined through equation Q below. 
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where 



M 12 (V avg ) = m Vlm - 25 - 51og(U avg /ffo 



Here, $(M) represents the differential galaxy luminosity 
function for the sample and Dq the projected separation 
(in Mpc) at some chosen fiducial redshift Vf- Mn m = 
Mi2(Vf) is a constant for a given Vf, and mi; m is the 
apparent-magnitude limit of the sample. 

This scaling of the linking length compensates for the 
bias that would otherwise be introduced due to the vari- 
ation in sampling of the luminosity function with red- 
shift. There is much debat e on how and whether or 
not to s cale V T (e.g., HG82. iNolthenius fc Wmtelll987t 
lFredenH lT9^5hM. If one assumes uniform density, sim- 
ple scaling arguments show that the velocity is simply 
proportional to the radius, suggesting that Vl should be 
scaled in the same manner as Dl- Such a scaling would 
include unwanted interlopers at large values of Vl and 
thus introduce an unwanted correlation between velocity 
dispersion and redshift. The density profiles of galaxy 
clusters, however, are usually better described by the 
isothermal-sphere approximation; in this case the veloc- 
ity dispersion is independent of the size of the cluster. 
It follows therefore that by setting Vl to a reasonable 
fixed value, we will minimize the number of interlopers, 
but not bias the algorithm against finding distant groups. 
Hereafter we set 

V L = V (2) 

The choice of Dq determines the minimum density con- 
trast of identified groups, which can be estimated using 
equation © below (HG82). 
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3. THE SAMPLE 



The first available sample of the 2MRS galaxy catalog 
contains positions, redshifts and magnitudes for 23090 
galaxies selected from the XSC. The targets were by se- 
lected by introducing a cut on the corrected magnitudes 
of objects in the XSC of K < 11.25 mag (the apparent 
magnitudes had previou sly been corrected fo r extinction 
using the dust maps of Schlcge l et al.lfl998)) . This cat- 
alog is complete, bar 40 galaxies, for galactic latitudes 
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|6| > 10° between galactic longitudes 330° and 30°, and 
\b\ > 5° for other longitudes. 

Below, we discuss a simple flow-field model applied to 
provide improved estimates of the distances to the galax- 
ies (see t!3.f \ , In ^ 13.21 we discuss the method used to 
populate the galactic plane with random galaxies to pre- 
vent any artifacts arising from the significantly reduced 
observed number density of galaxies behind the plane. 
We briefly discuss the assumed form of the luminosity 
function of the sample in tJ3.3l and, in tJ3.4l consider the 
completeness of the sample in redshift-space. 

3.1. Distance estimates 

Locally redshifts do not provide a reliable indication 
of distance because of distortions to the local velocity 
field due to infall onto concentrations of mass. Although 
the clustering algorithm is independent of the observer's 
frame of reference, it is essential to have reasonable esti- 
mates of the distances to the galaxies in order to compute 
the linking parameters, D\2 and Dl, as well to accurately 
estimate the luminosities of the galaxies. 

We apply t he b asic flow-field model described in 
iMould et alJ 1)200(1) to account for the local dis- 
tortions to the velocity field. This prescription 
first corrects the reference frame to the LG frame 
( Yahil . Tammann. fc Sand agdll97 7|. corroborated by th e 
more recent work of Courteau & van den Berghl 119991. 
then adjusts the redshift-inferred distances of galaxies 
near Virgo, Shapley and the GA region as follows: All 
galaxies within 12° of the center of Virgo with heliocen- 
tric redshifts less than 2500 km s _1 are placed at the red- 
shift of Virgo (plus a random velocity, drawn from a gaus- 
sian distribution with a standard deviation of 20 km s , 
to avoid artifacts in the group properties occurring from 
galaxies with identical redshifts). All galaxies within 
10° and 2000 km s _1 of the GA are placed at the red- 
shift of the GA (plus scatter) and all galaxies within 
12° and 3000 km s _1 of Shapley are placed at the red- 
shift of Shapley (plus scatter). The corrected velocities 
are then used in place of the heliocentric velocities when 
computing distances only. To infer the distances we as- 
sume Hubble's law is valid to the completeness limit of 
the 2MRS catalog, 3 using a Hubble constant Hq = lOO/i 
km s _1 Mpc _1 , where we assume h — 0.73 when a specific 
value is required. Th is value is chosen ba sed on the three- 
year WMAP results l|Spergel et all2006|) . h = 0.73±0.03. 
In the very local universe (i.e. where corrected distances 
are less than 3/i -1 Mpc) we give galaxies an indicative 
distance of 3/i -1 Mpc. The velocities used in computa- 
tion of Via, etc. are the heliocentric velocities reported 
in the 2MRS catalog. 

3.2. Filling in the galactic plane 

The 2MRS catalog is currently incomplete near the 
galactic plane (|6| < 10° between galactic longitudes 330° 
and 30°, and \b\ < 5° for other longitudes). With a signif- 
icantly reduced number density of galaxies observed be- 
hind the galactic plane, a structure that spans the plane 
will not be identified by the clustering algorithm. Sim- 
ilarly, structures that are visible in part above or below 

3 At the limiting redshifts of the galaxies analyzed, the difference 
between distances computed using a ACDM cosmology and simply 
assuming Hubble's law is less than 5%. 



the plane may not be identified as groups, and, even if 
they are, a bias will be introduced in the number density 
of groups with centers just above or below the plane. Any 
flow-field model derived from such a group catalog will 
suffer from these biasing effects; we therefore attempt to 
minimize these effects by randomly populating the sam- 
ple to enhance the galaxy number density behind the 
galactic plane to reflect that observed above and below 
it. 

We follow a method similar to that of I Yahil et alJ 
l)1991|) : this method produces similar results to the 
more involved W iener reconstruction 4 discussed in 
lLahav et alJ l|1994l) . We first divide the catalog into bins 
spanning 10° in galactic longitude and lO/i" 1 Mpc in dis- 
tance. For galactic longitudes ranging from 330° to 30° 
(masking the bulge) we now consider bins further bound 
by the lines |6| = 10°. Sampling from the adjacent bins 
(10° < \b\ < 20°), we populate the bulge with N galax- 
ies drawn at random from the galaxies in adjacent bins; 
these galaxies are placed at random latitudes and a nor- 
mal scatter of 20 km s _1 is introduced in the velocity to 
prevent artifacts in the group properties arising due to 
galaxies at identical redshifts. N is calculated by draw- 
ing a random normal deviate from a distribution with a 
mean equal to the number of galaxies in the two adjacent 
bins (above and below the plane), then subtracting the 
number of galaxies already present within the bin. For 
other galactic longitudes, the latitudes 5° < |6| < 15° are 
used to populate the bins with |6| < 5°. 5 The catalog, be- 
fore and after population of the galactic plane, is shown 
in Figure 2] The population of the plane generated an 
additional 2076 galaxies. 

3.3. Luminosity function 

The K-band luminosity function utilized in algorithm 
is parameterized in terms of a function of the form 
ISchechterlllT97ft. 

S(M) - 0.41n(10)$n0 a4 ( a+1 )( M *- M ) 
xexp[-10°- 4 ( M *- M )] 
We use the values reported in lHuchra et al.l ((2005b) , 
a = -1.02 

M* = -24.2 (4) 

$* = 1.08 x 10" 2 h 3 Mpc" 3 

which have been computed using the galaxies in the 
2MRS catalog with galactic latitudes, |6| > 10°. 

3.4. Completeness in redshift-space 

Due to the nature of flux limited surveys, the number 
density of galaxies observed at sufficiently high redshifts 
will tend toward zero. At these highest redshifts, the 
linking lengths used in the algorithm become so large 
that the majority of identified groups will likely be spu- 
rious. For the purposes of building a flow-field model, the 
groups at the highest redshifts will have the smallest af- 
fect on local dynamics, thus we choose to limit the group 

4 See lErdogdu et alJ §006), and Rassat et al. (2006), in prepa- 
ration for discussions on the Wiener reconstruction of the 2MRS 
sample. 

5 In this case, we set the mean of the normal distribution from 
which N is drawn to half the number of galaxies in the adjacent 
bins. 
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Fig. 1. — Galaxies in 2MRS catalog shown in a Mollweide projection in galactic coordinates. The top panel shows the catalog before 
the plane was populated. The bottom panel shows the catalog including the addition of the randomly-generated galaxies. The solid line 
indicates the region that was populated. (The color in the electronic edition indicates the measured redshift of the galaxy.) 
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Fig. 2. — Selection function of the 2MRS survey. The data points 
(shown as small dots) represent the cumulative number of galaxies 
as a function of (estimated) distance. The data are fit with a curve 
of the form of equation J2J using a least-squares fit (solid line). 
The dotted lines show the 5<r-errors from Poisson statistics. The 
data departs from the fit at the both smallest and largest distances 
shown on the plot. The dashed line shows the derivative of this 
curve (the selection function) in arbitrary units. The histogram 
contains the binned data shown in the same arbitrary units as the 
selection function for purposes of comparison only. The maximum 
and half-maximum values of the selection function are indicated. 

catalog to a redshift inside which the catalog is reason- 
ably complete. Figure El shows the cumulative number of 
galaxies, iV(< D), as a function of (estimated) distance, 
D. The points have been fitted with a curve of the form 

where the best fit parameters are No = 2.57 x 10 4 , 
a = 2.10, S = 106 Mpc, (3 = 0.881, b = 3.94. The 
derivative of equation (0) represents the selection func- 
tion of the survey, N(D), where the number of galaxies 
observed with estimated distances between D and D+dD 
is given by N(D)dD. We choose to cut the group catalog 
at the distance where the selection function falls to half 
its maximum, D cut = 140 Mpc, illustrated in Figure [3 
The entire 2MRS K < 11.25 dataset plus the galaxies 
generated in the population of the galactic plane will be 
used to create the group catalog, but the catalog will then 
be truncated, excluding groups with mean estimated dis- 
tances greater than D cut . 

4. PARAMETER CHOICES 

In this section we justify the choice of linking parame- 
ters used in the construction of our group catalog. Any 
group catalog produced from the remaining data sample 
will contain minimal biasing effects at the highest and 
lowest redshifts as well as minimal edge effects across 
the galactic plane. There remains, however, a choice of 
the parameters that specify the minimum density con- 
trast of detected groups. These parameters, Do and Vq 
from equations and J5J) above, must be chosen in a 
somewhat arbitrary fashion. 

There is no perfect choice of these parameters that will 
allow us to identify only groups which are gravitationally 
bound. In any choice we make, some bound systems may 



be divided and unbound galaxies will be present in some 
groups. For very large values of both D and Vb, the algo- 
rithm will associate all of the galaxies into a single group. 
Conversely, if we choose sufficiently small Dq or Vq, the 
algorithm will divide s ubstructures within real cluster s 
into multiple systems IfRamclla. Pisani. fc Oellerl ll997"l . 
eventually separating each galaxy into its own group. It 
is clear, therefore, that a suitable parameter choice will 
lie between these extreme cases. The method of choos- 
ing the specific values of the parameters must still re- 
main arbitrary; in order to be able to infer properties of 
the Universe (e.g. the matter density parameter) from 
the catalog it is unwise to calibrate the algorithm us- 
ing simulations based on a set of defined initial assump- 
tions as this would bias our results towards recovering 
these initial values. 6 It is obvious that there will be a 
choice of parameters that lie between these extreme val- 
ues which maximizes the number of groups produced. It 
is therefore reasonable to use a method of maximization 
to determine the choice of linking parameters, with no 
alternative method available that does not have similar 
shortcomings. 

At this point, we must consider the size of the group 
we choose to maximize. We choose to ignore binaries 
in our definition of groups as previous work has shown 
such systems identified using percolation algorithms to 
be un bound in the majority of cases (e.g lDiaferio et all 
1999). We consider the parameters obtained when max- 
imizing the number of groups of G or more members for 
3 < G < 20 as described below. We choose to set Vp = 
1000 km s _1 following HG82. Figure[3]shows the number 
of groups containing 3 or more members in Dq-Vq space. 
In Figure |2f a) we explore the parameter space on the 
intervals D = [0, 10] Mpc, V = [0, 2000] km s -1 . We 
then attempt to maximize the number of groups obtained 
by the following method: we divide the region spanning 
-> 10 Mpc in D and -> 2000 km s -1 in V into a 9x9 
grid and search for the combination of parameters that 
produces the largest number of groups. We then change 
the range of the Do and Vo parameters spanned to coin- 
cide with a 3x3 grid (as far as possible) centered on the 
values of Do and Vq that produced the largest number of 
groups. We divide this region into a 9x9 grid and itera- 
tively repeat the procedure until the desired accuracy of 
the parameters is reached. This procedure is illustrated 
in Figure I^Jb). We repeat this maximization procedure 
for all values of G between 3 and 20. The obtained value 
of Do rises gently with G, while the scatter in obtained 
values of Vq increases rapidly with G. As observed in 
Figure Et a ); the number of groups depends highly on 
Do, but is comparatively insensitive to Vq. Since the 
velocity dispersion of a cluster is not expected to de- 
pend on the size of the cluster (see the discussion in 
above) it is not sensible to consider large values of Vq (i.e. 
> 1000 km s _1 ) as this will introduce members that are 
not gravitationally bound and thus propagate errors into 
the mass estimates. The fraction of interlopers present 
in a group will also increase with both Do and Vbj thus 
the most sensible choice of parameters to minimize inter- 
lopers and reduce the scatter in Vb corresponds to G = 3. 

The maximum number of groups of 3 or more members 

6 See Crook et al. 2007a, in preparation, for further investiga- 
tion. 
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Fig. 3. — The number of groups of 3 or more galaxies obtained 
algorithm has been executed for each pair of parameters on a 25 X 
maximization routine discussed in the Jl]of the text. 

(1538) is obtained for the values (Do, Vq) = (1.63 ± 
0.03 Mpc, 399 ± 8 km s _1 ), corresponding to the density 
contrast Sp/p = 12. 

In an analysis of the northern CfA redshift survey 7 
(hereafter CfANI. IRamella et al] l)1997f) show that the 
group properties are statistically stable for values of 
density contrasts Sp/p > 80, where they scale Vl in 
a similar manner to Dl and choose Vq = 350 km s _1 . 
iDiaferio et al.l l)1999j) apply a similar choice of parameters 
to mock CfA surveys based on A-body simulations and 
conclude that 80% of groups with 4 or more members are 
true virialized systems, wher eas 40% of triplets are not, 
confir ming the hypothesis of lRamella. Celler. fc Huchral 
(1989). As Figure O^a) shows minimal variation in the 
number of groups produced wi th changing Vn comp ared 
to changing Do, the findings of lRamella et alJ l)1997(l are 
applicable to this study. We will proceed to analyze the 
groups produced at both the values of D and Vq that 
maximimize the number of groups of 3 or more members , 
as well as the values suggested bv IRamella et alJ l|1997Tl 
(i.e. Do = 0.89 Mpc, which corresponds to Sp/p = 80, 8 
and Vq = 350 km s" 1 ). 

5. GROUPS 

We present the results of applying the group- finding al- 
gorithm (SJ2J to the 2MRS catalog subset (fJ3J using both 
pairs of parameters discussed in above. The group 

7 This is a subset o f the extended CfA r e dshift survey 
Ide Lapparent et all Il99ll IGeller fc Huchral ITflSflt IHuchra et al] 
119901119951) 1 

8 The smallest allowed density contrast is chosen to minimize 
the probability of splitting the richest systems. 
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(b) 

as a function of the parameters Do and Vq. In (a), the clustering 
grid, (b) contains a graphical representation of the execution of the 

catalogs are presented in Tables lA"TTlA6l in the appendix 
(see the electronic edition for the complete catalogs) . We 
provide an overview of the catalogs in i J5.ll below, then 
discuss the identified groups and contrast the two cata- 
logs in i)5.2l We present the properties of the obtained 
groups f H5.3J) and discuss the reliability of the clustering 
algorithm f H5.4(l . 

5.1. Overview 

The catalog produced using the parameters (Dp , Vq) — 
(1.63 Mpc, 399 km s _1 ) is presented in Table lAll These 
parameters produced the maximum number of groups 
of 3 or more galaxies, and correspond to a density con- 
trast Sp/p = 12; this catalog will hereafter be referred to 
as the low-density-contrast (LDC) catalog. The catalog 
produced using the parameters (-D p, V q) = (0.89 Mpc, 
350 km s _1 ) is presented in Table lA"2l These parame- 
ters correspond to the density contrast Sp/p = 80; this 
catalog will hereafter be referred to as the high-density- 
contrast (HDC) catalog. 

Figure 01 shows the positions and sizes of all groups in 
the two catalogs. Figures Ufa) and^b) show the groups 
in equatorial coordinates; FiguresEfc) and^Jd) show the 
groups in galactic coordinates (available in the electronic 
edition only). The maps are shown in Mollweide pro- 
jection, which preserves the area of the structures on the 
surface of a sphere but distorts their shape, especially 
close to the poles. The plots show ellipses that have 
the properties of the groups discussed in t |5.3.2l below, 
and are only representative of the shape and size of the 
group. The Local SuperCluster (LSC) has been clearly 
identified in Figure EJa) as the large structure in the 
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Fig. 4. — Groups (scaled by angular size) identified by clustering algorithm in 2MRS catalog. Figures (a) and (b) are shown using a 
Mollweide projection in equatorial coordinates, centered at 12 ,i 00 m . Figures (c) and (d) (available only in the electronic edition) are shown 
using a Mollweide projection in galactic coordinates. Each group is plotted as an ellipse wi th the major axis, position angle and axis-ratio 
of the group. The ellipses are transformed from the x-y coordinate system discussed in H5.3.2l to the appropriate map projection. The 
color of the ellipse represents the group's mean redshift. The galactic plane is shown by the dotted line in Figures (a) and (b). Figures (a) 
and (c) show the groups in the LDC catalog (Sp/p = 12); the LSC has been clearly identified as the central structure. Figures (b) and (d) 
show the groups in the HDC catalog (Sp/p = 80); we observe the effect of increasing the minimum density contrast required to identify 
groups: the LSC has been broken into several constituents. The distortions created by the map projection are maximized near the poles, 
and enhanced for larger structures. The seemingly strange shape of the LSC in Figure (c) is the result of mapping an ellipse in the x-y 
coordinate system onto this projection, and is only partly representative of the true shape of the structure. 
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center of the figure. When applying the algorithm with 
the higher minimum density contrast, this structure is 
split into several constituents as shown in Figure Elb). 
(The same result is evident in Figures HJc) and Eld) in 
the electronic edition, however the LSC encompasses the 
pole of the coordinate system in this case. The distor- 
tions due to the map projection are therefore enhanced 
in these plots, and the shape is less representative of the 
true shape of the LSC. The area occupied by the LSC 
in Figure 0Jc) is the same as that in Figure Eld) . The 
constituents that have been merged to form the LSC are 
clearly visible in Figure Eld).) 

There is an apparent enhancement in the number of 
groups with large angular sizes near the galactic plane 
(see Figures Elb) and Eld), the latter is available in the 
electronic edition only). There are 5 groups shown as 
ellipses with major axes greater than 5° and with cen- 
ters inside |6| < 10°. Of these, 2 have only 3 members (of 
which 2 are genuine galaxies from 2MASS XSC) and only 
1 out of the remaining 3 groups has more than 10% of its 
members randomly generated in the population of the 
plane; the remaining 2 are the only groups with more 
than 5 members, thus the large apparent sizes of all 5 
groups in the figure is due to their proximity. We con- 
clude that the observed enhancement is therefore not an 
artifact of the population of the plane. 

Figure shows the same groups as Figure EJ however 
the area of each ellipse is proportional to the number 
of members in the group, rather than the square of the 
group's angular size. The areas are normalized such 
that the major-axis (see H5.3.2|) of the largest group ap- 
pears as 75% of its true angular size. These figures are 
no longer dominated by the foreground groups that have 
the largest angular sizes, but show how the groups are 
distributed across the sky out to the redshift limit of the 
sample. 

5.2. Identification and overlap 

Due to the the nature of the group finding algorithm, 
we expect to find that all of the galaxies assigned to 
groups in the HDC catalog will also be assigned to groups 
in the LDC catalog, however the converse is not nec- 
essarily true. We consider all the galaxies assigned to 
groups in the HDC catalog and determine the corre- 
spondence between groups in the two catalogs. The six 
largest groups in the LDC catalog are plotted in Fig- 
ure [HJ the corresponding groups in the HDC catalog 
are also shown. The largest group in the LDC cata- 
log contains 810 members (of which 2 were randomly 
generated in the population of the galactic plane). This 
group is the result of merging several smaller groups in- 
cluding Virgo, NGC3607, NGC4105, IC764, NGC5746, 
NGC3190, NGC5846, and NGC4038 Clusters, as well as 
the M81 group, and corresponds to the LSC. In the HDC 
catalog, most of these groups have been identified indi- 
vidually; in fact Virgo has been split into two groups, 
containing 298 galaxies and 123 galaxies, respectively. 
Eridanus, Fornax I, Dorado and NGC2280, NGC1433, 
NGC2559 Clusters merge to form the second largest 
group in the LDC catalog (302 members, including 22 
simulated) ; again these were identified individually in the 
HDC catalog. A426 and A347 of the Perseus-Pisces su- 
percluster make up the third largest group in the LDC 
catalog, containing 301 galaxies, but these were identi- 



fied as two separate groups in the HDC catalog. Hydra 
(A1060) is identified as the fourth largest with 241 mem- 
bers, and Norma (A3627, the GA) is identified as the 
fifth with 217 members (42 of which were randomly gen- 
erated). Centaurus (A3526) was identified as the sixth 
largest group (202 members). This exercise demonstrates 
that the correct choice of parameters used in the group 
algorithm is entirely dependent on the size of the struc- 
tures that are desired. 

The remaining identifications and correspondence be- 
twe en th e catalogs are shown in the appendix ( Tables lA"3l 
and lA4(l . where we consider only groups cont ainin g 2 5 or 
more members in the HDC catalog. Tables lA~5l and lA~6l 
contrast the group assignments of individual galaxies be- 
tween the two catalogs. It may be surprising that the pa- 
rameters chosen to maximize the total number of groups 
actually merge several of the large groups, hence appar- 
ently reducing the total number of groups. Although not 
obvious, this result is not unexpected because the larger 
linking length will allow many smaller groups to be iden- 
tified that do not exist in the HDC catalog. The latter 
association is generating more groups than are removed 
by the merging of largest groups. It is likely that the 
groups in the LDC catalog contain a higher fraction of 
interlopers than the groups in the HDC catalog, however 
it is evident that the LDC catalog identifies the largest 
structures on the sky. This suggests that the LDC cata- 
log will be the better candidate for the basis of a flow-field 
model, as some of the largest structures are fragmented 
in the HDC catalog. This is explored further in follow-up 



5.3. Group properties 

In this section, we discuss the properties of the LDC 
and HDC catalogs. Estimates of the velocity dispersion, 
size, mass and luminosity are discussed in S15.3.1I we ob- 
tain estimates of the axis-ratio and position angle of the 
groups in H5. 3. 21 below. 

5.3.1. Basic properties 

The properties of the LDC and HDC catalogs are sum- 
marized in Tabled in further analysis, we only consider 
groups with 5 or more members that are also present 
in the 2MRS catalog (referred to as genuine hereafter), 
as opposed to those generated in the population of the 
galactic plane. We do this in an attempt to exclude 
groups with a high fraction of interlopers in our analysis. 
We provide two estimates of the mass of the groups. We 
first compute the virial mass of the group, My, 



M v = 



3ir J2p Rp 
T G 



(6) 



where G is Newton's constant, op is the projected veloc- 
ity dispersion, 



2 ZW-Vg) 2 

CTp= JV-1 

and Rp is the projected virial radius, 

N(N - 1) 



R P = 



i>3 ij 



(7) 



(8) 



9 See Crook et al. (2007b), in preparation. 
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Fig. 5. — Groups (scaled by number of members) identified by clustering algorithm in 2MRS catalog shown in equatorial coordinates. 
Figures (a) and (b) are shown using a Mollweide projection in equatorial coordinates, centered at 12 00 m . Figures (c) and (d) (available 
only in the electronic edition) are shown using a Mollweide projection in galactic coordinates. Figures (a) and (c) show the groups in the 
LDC catalog (Sp/p = 12). Figures (b) and (d) show the groups in the HDC catalog (Sp/p = 80). Each group is plotted as an ellipse with 
the position angle and axis-ratio of the group. The areas of the ellipses are proportional to the number of members in the group, scaled 
suc h that the major-axis of the largest group is 75% of its true size. The ellipses are transformed from the x-y coordinate system discussed 
in i|5.3.2l to the appropriate map projection. The color of the ellipse represents the group's mean redshift. The galactic plane is shown by 
the dotted line in Figures (a) and (b). 
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the right-hand figure at (13h, +10°) is partially masked by the high density of points. Middle: Eridanus, Fornax I, Dorado, NGC2280, 
NGC1433, NGC2559 Clusters (LDC Group #391). Bottom: Perseus-Pisces (A426, A347) (LDC Group #229). 2nd page: Top: Hydra 
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TABLE l 

General properties of the group catalogs 



Property 


LDC Catalog 


HDC Catalog 


Sp/p 


12 


80 


DO (Mpc) 


1.63 


0.89 


VO (km s" 1 ) 


399 


350 


No. of Singles 


5548 (27.2%) 


9608 (47.2%) 


No. of Binaries 


1397 (13.7%) 


1710 (16.8%) 


No. of Groups of 3+ 


1538 (59.1%) 


1258 (36.1%) 


No. of Groups of 10+ 


203 (30.5%) 


113 (13.3%) 


No. of Groups of 50+ 


17 (10.7%) 


8 (4.7%) 


Mean no. per group 


7.84 


5.84 


Standard Deviation 


26.25 


11.85 


Min/Max per group 


3 / 810 


3 / 298 



Note. — The values in parentheses represent the percentages of galaxies that fall into this category. 

Vg is the mean group velocity, and V% is the line-ol-sight 
velocity of the i th member. N is the total number of 
galaxies in the group and Rij is the projected separation 
between two galaxies, defined in terms of their angular 
separation, 6^, through 



% __tan(- 



Due to the biases in the virial mass estimator (e.g. 
iBahcall fc Tremainel 11981^) . we also calculate the pro- 
jected mass estim ator, Mp ijBahcall fc Tremainel Il98lt 
IHeisler et aUllflS^. 



M P = 



/PM 



nG(N - 7) 



(9) 



where Sj is the offset of the i member (in p hysical units) 
from the center of the group. Following Hcislcr et al. 
(1985f), we set 7 = 1.5 and / PM = 10.2. 

We estimate the total isophotal K-band luminosity, Lk 
of each cluster from the observed luminosity using equa- 
tion ijTUj) below. 



J K 



1 - 



7 (a + 2, L lim /L*) 



L 



obs 



(10) 



T(a + 2) 

where L b s is the total isophotal observed K-band lumi- 
nosity and L]x m is the limiting observable luminosity at 
the distance of the cluster, d c , 

L Vlm = io o - 4 ( M 0. K - m ii'»+ 25 + 51o s(' it =/ M P c ))L Q (11) 

a and L* take the values used quoted in equation (0J 
above, given 

L* = io°- 41 °s( m o.k- m *) l 
and 7(777, x) is the lower incomplete gamma function, 



7(m, x) = 



dt 



and 



F(m) = 7(771, 00) 



We set the K-b and magnitude zero point, Mq,k = 3.39 
ij.Tohnsonl I1966D . Note that we have not applied an 
isophotal correction to the luminosities, thus the lumi- 
nosities presented in this paper are lower than the total 
K-band luminosities. 

Table |2 contains the median properties of the groups 
in the catalog that have at least 5 genuine members. 



Figure [7\ shows the fraction of groups, containing at 
least 5 genuine members, as a function of velocity dis- 
persion, projected virial radius, mass and mass-to-light 
ratio (using both virial and projected mass estimates). 
As expected, the distribution of virial radii widens and 
increases with the larger choice of Dq. Similarly, the ve- 
locity dispersions increase with the larger choice of Dq 
and Vq due principally to the change in Vq ; re-calculating 
the groups for various values of Dq while keep Vq constant 
has little impact on the distribution of velocity disper- 
sions. The net effect is expected to increase the estimated 
virial mass of the groups, as observed. 

We use the median mass-to- light ratio to obtain a value 
of f2 m , the ratio of the matter-density, p m (assuming 
all the mass is contained within galaxy clusters) to the 
critical density, p Cl 



Pc = 



ml 

8nG 



(12) 



We use the K-band luminosity density, Ck, 

C K = $*i*r(a + 2) 

where a, L* and $* take the values given in equation 
in conjunction with the median mass-to-light ratios 
to estimate p m , 



,Me 



K 



(13) 



where E refers to the method of mass estimation (i.e. 
Virial or Projected). For the luminosity function and 
magnitude zero point used in this paper, the critical 
density corresponds to ( Mr/Lk) = 2>22M ( 7 - 1 /Lp t . Sim- 
ilar result s were found bv iKochanek et alTl)2001|) and 
iBell et aTl l)2003l) . adjusting for the authors' isophotal 
corrections where appropriate. 
The obtained values are given in Tabled Spergel et alJ 



2006) obtain Vl r 



0.238. 



-0.013 
-0.024 



assuming h = 0.73; this 



value only agrees, at the la-level, with the value we ob- 
tained using the projected mass estimator in the LDC 
catalog (Sl m = 0.2291qq^2i 1°" errors). The virial mass 
estimates predict f2 m to be too small. For this reason, 
we will use the projected mass estimator as opposed to 
the virial mass estimator in further analysis. The fact 
that the HDC catalog predicts a value of VL m significantly 
smaller than that obtained using WMAP suggests that 
we are missing a significant fraction of the mass of the 
cluster in our estimate; selecting groups based on the 
density contrast Sp/ p = 80 is causing us to underestimate 
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Fig. 7. — Properties of groups with at least 5 genuine members. We contrast the properties of the groups in the LDC catalog (solid 
line) with those in the HDC catalog (dashed line). Top left: variation in line-of-sight velocity dispersion across the groups. Top right: 
shows the projected virial radius of each of the groups. Middle: Virial (left) and projected (right) mass estimates of the groups. Bottom: 
Mass-to-light ratios of the groups using the virial (left) and projected (right) mass estimators. We assume h = 0.73. 
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TABLE 2 

Median properties of groups with 5 or more genuine members 



Property LDC Catalog HDC Catalog 

<rp (km s" 1 ) 197 (183, 206) 183 (166, 193) 
Rpv (Mpc) 1.71 (1.58, 1.85) 0.97 (0.90, 1.04) 
\og[M v /M Q ] 13.79 (13.72, 13.90) 13.54 (13.46, 13.60) 
log[M P /M ] 14.05 (13.98, 14.10) 13.66 (13.57, 13.75) 
log[(M v /L K )/(M Q /L Q )] 1.66 (1.59, 1.72) 1.49 (1.41, 1.55) 
log[(M P /L K )/(M e /L e )] 1.86 (1.81, 1.93) 1.63 (1.56, 1.69) 
Q M ,v 0.14 (0.12, 0.16) 0.10 (0.08, 0.11) 
Q M ,p 0.23 (0.20, 0.27) 0.13 (0.11, 0.16) 

Note. — The median values are shown, with 99% confidence levels in parentheses. We compute the confidence levels by drawing an 
equally-sized sample from the observed distributions and computing the median values 5,000 times. 
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TABLE 3 

Velocity-dispersion- and M / L— distance relations 



50 100 150 

Position Angle (degrees) 



Fig. 8. — Position angles of groups with 5 or more genuine 
members. The figure shows the fraction of groups with position 
angles in the specified intervals (binned by 10°). The groups in 
the LDC catalog are represented by the solid line, and the HDC 
catalog by the dashed line. 

the median mass-to-light ratio of groups. This suggests 
that the dark matter halos extend beyond the Sp/p — 80 
density-contrast contour that was inferred from luminous 
matter. 

5.3.2. Orientation and ellipticity 

We include in the catalog a measure of the axis-ratio, 
position angle and semi-major axis of the groups con- 
taining 5 or more genuine members, calculated using the 
following method: We rotate the coordinates such that 
center of the group 10 lies along the z-axis. We measure 
the angle of each galaxy from the z-axis, O, as well as its 
azimuthal angle ip, then define 

z = ecos(V>) (14) 
y = esin(V>) (15) 

such that the positive y-axis points north and the posi- 
tive x-axis points east. 11 This definition is chosen such 
that the shape of the group is not distorted under the 
projection onto a plane. 

We rotate the axes to some angle, <fi, and define the 
coordinates of the galaxies in the rotated frame as (ii, 
yi). We choose the value of </> that minimizes S^y 2 . We 
compute the 75 th -percentile values of \x\ and \y\ (275 and 
j/75) and record their ratio, ij, (0 < r\ < 1) as a measure 
of the axis-ratio of the group. We also record the larger 
of X75 and j/75 as a measure of the semi-major axis of the 
cluster, a, as well as the angle of rotation of the semi- 
major axis from north toward east (the position angle of 
the group, 0). We verify that these angles are approxi- 
mately uniformly distributed by showing the number of 
galaxies as a function of position angle in Figure |H1 

The above properties of the groups are illustrated 
graphically, by using ellipses with the same semi-major 

10 We define the center as the mean position of the galaxies in 
cartesian coordinates, assuming the galaxies lie on the surface of a 
unit sphere. 

11 Note that the mean values of x and y are not strictly zero 
but, in the analyzed data, are sufficiently small that an iterative 
centering procedure is not required. 



Parameter 


LDC Catalog 


HDC Catalog 


a (10-* Mpc" 1 ) 


1.0 ±0.3 


1.2 ± 0.3 


(3 


2.21 ±0.03 


2.17 ±0.03 


p (10-* Mpc" 1 ) 


-7.0 ±0.6 


-5.8 ±0.6 


q 


2.43 ± 0.05 


2.02 ± 0.05 



Note. — Values correspond to the parameters in equations H16I 
and 1171 that minimize their respective x 2 ~statistic. 

axes, axis-ratios and position angles in FigureQ]and with 
the same axis-ratios and position angles in Figure [S] For 
groups with 3-4 members, a circle is drawn with an an- 
gular radius equal to the 75 th -percentile mean offset. 12 

In Figure IHI we show the galaxies that are associated 
by the clustering algorithm to form the six largest groups 
in the LDC catalog. The corresponding ellipses have 
been overlayed on this plot. 13 The figure also shows the 
corresponding groups identified in the HDC catalog. It is 
evident that the higher density contrast used in the latter 
choice of parameters splits the large structures identified 
when choosing a lower-density contrast. 

5.4. Reliability of the algorithm 

In this section we discuss the verifications performed 
to ensure that the groups obtained are consistent with 
both expectation and the literature. In 15.4.11 we exam- 
ine the distance-dependence of the velocity dispersions 
and mass-to-light ratios of the groups. We compute the 
mass functions of the group catalogs and compare them 
with expectation in l5.4.2l Finally, we compare the 2MRS 
group catalogs di rectly with the UZC-SSRS2 and CfAN 
group catalogs in 15.4.31 

5.4.1. Variation with distance 

Figure [5] shows the velocity dispersion, erp, of the 
groups as a function of distance. We fit a curve of the 
form 



a P = 10( aD +V km s" 



(16) 



to the data, where the best-fit parameters are given in 
Table|21 The large scatter and very small correlation ob- 
served in Figure [5] (the mean velocity dispersion changes 
by a factor of 0.635- between [20, 40] Mpc and [120, 140] 
Mpc in the HDC catalog, and 0.35tr in the LDC cat- 
alog 14 ) demonstrates that there is minimal bias intro- 
duced in the velocity dispersion of groups with distance; 
this was desired in the construction of the algorithm (see 
0. Had we chosen to scale Vo with distance, we would 
expect the velocity dispersions of the most distant groups 
to be larger than observed in this figure. Since the corre- 
lation is already slightly positive, scaling Vo would have 
introduced a more significant bias with distance. 

Figure^|shows the mass-to-light ratios (computed us- 
ing the projected mass estimators) as a function of dis- 
tance for both pairs of parameters. We fit a curve of the 

* 2 T his va lue is also reported under the column titled a in Tables 
lATlandlA^l 

13 Note that the shapes of the ellipses have been distorted due 
to the choice of coordinate system. 

14 <t represents the average standard deviation, weighting the 
standard deviation of the velocity dispersions in each interval 
equally. 
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Fig. 9. — Velocity dispersion of groups as a function of distance. The left panel shows the groups in the LDC catalog (Sp/p = 12). The 
right panel shows the groups in the HDC catalog (5p/p = 80). Only groups with 5 or more genuine members are included. The solid line 
shows a linear fit to the data as described in the text. 
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Fig. 10. — Mass-to-light ratio of groups as a function of distance. The plots show the ratio of the projected mass estimate to the corrected 
K-band luminosity of the cluster. The left panel shows the groups in the LDC catalog (Sp/p = 12). The right panel shows the groups in 
the HDC catalog (Sp/p = 80). Only groups with 5 or more genuine members are included. The solid line shows a linear fit to the data as 
described in the text. 



form 

Ml = i^ipD+q) M®. (17) 

Lk L© 

to the data, where the best-fit parameters are given in 
Table |3] The lower limit to the mass-to-light ratios com- 
puted as a function of distance remains approximately 
constant, while the upper limit decreases with distance, 
giving rise to the negative slope. Due to the nature of the 
flux-limited sample, at the largest distances we are not 
sensitive to (intrinsically) faint objects, therefore we are 
preferentially selecting the brightest groups. The scaling 
of the linking length, Dl, is designed to produce groups 
with a similar number of members at all distances. Since 
we are further correcting the luminosity of the groups 
to account for those galaxies to which the survey was 
not sensitive, we expect the mean luminosity of groups 
to increase with distance. As we have already shown 
that the velocity dispersions of the groups we find (and 
hence the estimated masses) are comparatively uncor- 



rected with distance, we would expect that we should 
miss those groups with high mass-to-light ratios at the 
largest distances, as indeed we observe in Figure EH To 
correct for this effect, one may introduce a scaling in the 
linking length in velocity space, Vl, with distance. This 
would increase the estimated group mass with distance, 
however such mass estimates would be based on groups 
containing many interlopers and thus not accurately rep- 
resent the mass of the group. Such shortcomings of the 
percolation algorithm will be discussed in more detail in 
follow-up work. 15 

5.4.2. Mass functions 

The large number of groups in the sample allows us 
to obtain an accurate estimate of the mass function for 
groups in the LDC and HDC catalogs. In this section we 
only consider groups with at least 5 genuine members at 
a distance of 10 /h Mpc or greater. We compute the mass 

15 See Crook et al. (2006a), in preperation. 
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TABLE 4 

Best-pit values for Press-Sohechter and Sheth-Tormen 
mass functions 



Parameter 


LDC Catalog 


HDC Catalog 


PS74 form 11 


p (10 1U M Q Mpc" 
log(M /M Q ) 


■ s ) 2.8 ±1.4 
14.50 ±0.13 


0.9 ±0.6 
13.93 ±0.18 


ST99 form b 


p (10 1U M Q Mpc" 
log(M /M Q ) 


J ) 3.6 ±1.9 
14.43 ±0.16 


1.3 ±0.8 
13.87 ±0.21 



Note. — Values correspond to the parameters in equation 1191 
that minimize the x 2 ~statistic. 
IPress fc Schechterl 119741) 

teheth fc Tormenl ll 99Hh . . 

function using the 1 / VVax procedure (e.g. lMartmez et al.1 
2002), whereby each group is weighted by the inverse 
of the maximum comoving volume, V mayL {Li), in which 
the group remains observable given the flux limit of the 
survey. Li is the luminosity of the fifth brightest member 
of the group. The differential mass function can be then 
computing as 



n(M) = 



\AI, 



-M\<AM 



(18) 



where Mi are the group masses and AM is the (variable) 
bin width. The results are shown in Figure ITT1 

We consider analyti cal differential mass fun ctions of 
the form suggested bv lSheth fc Tormenl <|1999ft (also see 
I Jenkins et al.ll2001^ . 



7/2 



x exp 



(2a) 



( M 

\M 



( M 

\M 



(19) 



where 7=1 
parameters a 



n/3) and we set h = 1. The choice of 
1, p = 0, A = 0.5 corresponds to the 
analytical pr ediction of | Press fc Schechterl fl97l (here- 
after, PS74). iSheth fc Tormenl l)1999|) suggest the alter- 
native choice parameters A = 0.3222, a = 0.707, p = 0.3 
(hereafter, ST99), which provide good agr eement with a 
subset of ./V-body simulations analyzed bv lJenkins et all 
(20Q3). 

We fit functions of the form of equation (|19fl to groups 
in the LDC and HDC catalogs using the parameter 
choices of both PS74 and ST99 (see Figure [H] and Ta- 
bic . The analytical descriptions are both good ap- 
proximati ons to the data ; simil ar conclusions were also 
drawn bv lMartmez et all l)2002|) . however we find the fit 
to the PS74 form produces a slightly smaller x 2_ statistic 
in both cases. 

We compare the ratio of the best-fit values for Mq in 
the PS74 form of the mass function to the value predicted 
using the simple arguments of PS74. Mo scales with the 
minimum density contrast according to 



9 



-2/7 



(20) 



and thus we expect the ratio of the determined values of 
M to be given by 

lQ g : ,mnr> = lo § 



M. 



(HDC) 




(<Wp)ldc 
(<Wp)hdc 



-2/7 



log (i)= 1.24 



The obtained values of Mq give log-ratio of 0.57 ± 0.22, 
providing agreement only at the 3cr-level with Press- 
Schechter theory. The two methods cannot be ex- 
pected to be within perfect agreement as the group- 
identification algorithm will find different groups as Vq 
is varied but Do is held constant. The computed esti- 
mates of the group masses will therefore vary, while the 
ratio predicted by the Press-Schechter treatment (which 
is sensitive only to the change in density contrast, thus 
Do) does not. 

5.4.3. Comparison with other group catalogs 

To verify the validity of the produced group catalogs, 
we compare the 2M ASS group catalogs w ith the UZC- 
SSRS2 group c atalog llRamella et alJ2002|) and the CfAN 
group catalog l|R.amella*'et*ldTlT997fl . The former is con- 
structed from partial versions o f the Updated Zwicky 
Catalog ( UZC: iFalco et al.lll999ft and the Sou thern Sky 
Redshift Survey fSSRS2: Ida Costa et al.lll998|) and cov- 
ers 37% of the sky; the latter covers 10%. The selec- 
tion criteria for the three catalogs result in completeness 
limits at different velocities. To compute the complete- 
ness in redshift-space, we use the same technique as M3.4I 
above, however in this case, we determine the velocity 
corresponding to the peak in the derivative of equation 
JHJ). When comparing the 2MASS group catalog with the 
other two, we cut both catalogs at the smallest of the two 
corresponding velocity limits (2MRS at 5697 km s _1 in 
the LDC catalog and 5350 km s _1 in the HDC catalog, 
UZC at 7115 km s" 1 , CfAN at 9390 km s" 1 ). 

For each of the groups present in the UZC-SSRS2 and 
CfAN catalogs, we search the 2MRS group catalogs for a 
group within a radius set by twice the sum of the virial 
radii of the group in 2MRS and the group in the compar- 
ison catalog, with mean velocities that differ by less than 
30%. We find that 86% of the groups in UZC-SSRS2 and 
76% of the groups in CfAN are present in the 2MRS LDC 
group catalog, whereas 78% of the groups in UZC-SSRS2 
and 69% of the groups in CfAN are present in the 2MRS 
HDC group catalog. There is reasonable agreement in 
both cases given the differing selection criteria employed 
in the three samples. 

We also compare the distributions of velocity disper- 
sion using the two-s ided Kolmogorov-Smirnov (K-S) test 
l)Chakravarti et all 11967ft . The cumulative fractions of 
the LDC and HDC catalogs are compared with the cu- 
mulative fractions of the UZC-SSRS2 and CfAN cata- 
logs (see Figure El- For comparison, we also show the 
velocity dispersions of groups in Sloan Digital Sky Sur- 
vey D ata Release 3 (SDSSDR3, iMerchan fc Zandivared 
2005) . The results of the test are shown in Table 
The velocity dispersions in both the LDC and HDC cat- 
alogs are consistent with the CfAN sample. However only 
the LDC catalog produces groups with velocity disper- 
sions consistent with those in the UZC-SSRS2 catalog, 
although the discrepancy with the HDC catalog is small. 
Any discrepancy with the CfAN catalog is not apprecia- 
ble due to the smaller number of groups in the compared 
sample of CfAN galaxies. 16 Neither the LDC or HDC 
catalogs are consistent with the SDSSDR3 group cata- 

16 This is accounted for in the Kolmogorov-Smirnov test. 
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Fig. 11. — Mass functions of the groups. The left panel shows the differential mass function of the groups in the LDC catalog (Sp/p 
= 12); the right panel shows that of the groups in the HDC catalog (Sp/p = 80). The masses are estimated using the projected mass 
estimator. 2cr error-bars are shown on the plots. Only groups at a distance of at least 10 Mpc containing at least 5 genuine members are 
included. The points have been fitted with a Press-Schechter (PS74) mass function (solid line) and a Sheth-Tormen (ST99) mass function 
(dotted line). The former fit produces the smaller x 2 -statistic in both cases. We assume h = 0.73. 
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Fig. 12. — Cumulative fraction of velocity dispersions of groups. The velocity dispersions in the UZC-SSRS2, CfAN and SDSSDR3 
group catalogs catalogs are compared with the 2MRS LDC catalog (left panel) and the 2MRS HDC catalog (right panel). The vertical 
lines indicate the corresponding D-statistic used in the K-S test. 



log. The discrepancy is expected due to the difference in 
the algorithms used to identify groups. 

The difference in the distributions of velocity disper- 
sion is due in part to the different scalings of the velocity- 
linking parameter. When scaling the linking length, the 
more distant groups are likely to have higher velocity 
dispersion. These groups may not be present in a group 
catalog derived by setting the parameter to a constant. 
The velocity dispersions of groups in the USZ-SSRS2 and 
CfAN catalogs tail off at -400 km s" 1 . The LDC cata- 
log is in reasonable agreement with this since the veloc- 
ity linking parameter is —400 km s _1 . The HDC catalog, 
however, contains a smaller fraction of groups with veloc- 
ity dispersions between 350-400 km s _1 than are present 
in the other catalogs. This is consistent with the above 
hypothesis since, in this case, the velocity linking param- 
eter is set to 350 km s _1 in this case. Without knowledge 



of the complete phase-space positions of the galaxies, we 
cannot determine whether the majority of groups with 
velocity dispersions greater than 350 km s _1 in the LDC, 
UZC-SSRS2 and CfAN catalogs are bound or not. 

We further examine the number density of groups 
as a function of redshift (see Figure I13|) . For com- 
parison, we include the number density of groups in 
SDSSDR3 and the 2dF GRS Percolation- Inferred Galaxy 
Group (2PIGG) catalog ijEke et al.l2004[) . although these 
surveys have very different selection biases and fainter 
flux limits, which result in a higher number density of 
groups than that detected in 2MRS. The 2MRS group 
catalogs contain an approximately constant number den- 
sity of groups to —80 Mpc, before following a power-law 
decay. The LDC catalog has slightly more uniform cov- 
erage out to this distance than the HDC catalog. The 
observed shape of the function is consistent between the 
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TABLE 5 

Comparison between 2MRS, UZC-SSRS2. CfAN and SDSSDR3 groups 



Property LDC Catalog 


HbC Catalog 


No. of 2MRS Groups 736 


630 


Comparison with UZC-SSRS2 


No. of UZC-SSRS2 Groups 444 


399 


P>-statistic a 0.070 


0.102 


P-valuc b 0.128 


0.011 


Comparison with CtAN 


No. of CfAN Groups 119 


107 


D-statistic a 0.128 


0.062 


P-value b 0.062 


0.862 


Comparison with SDSSDR3 


No. of SDSSDR3 Groups 204 


177 


D-statistic a 0.227 


0.122 


P-valuc b 10" 7 


0.029 



Note. — The table presents the results of two-sided Kolmogorov-Smirnov test on the comparison between the distributions of velocity 
dispersions of the groups in 2MRS catalog with the UZC-SSRS2, CfAN and SDSSDR3 catalogs. 

a D-statistic used in computation of P-value in Kolmogorov- 
Smirnov test. 

b P-value represents the probability that such a difference would 
be observed under the assumption that the two samples were drawn 
from the same parent distribution. We consider values of P < 0.05 
to indicate that the two samples were drawn from significantly 
different parent distributions. 
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Fig. 13. — Number density of groups as a function of distance. We show the number density for the LDC (diamond) and HDC (cross) 
catalogs, as well as the UZC-SSRS2 (triangle), CfAN (square), SDSSDR3 (open circle) and 2PIGG (filled circle) catalogs. (The symbols 
are colored black, blue, green, red, orange and yellow, respectively, in the electronic edition.) We assume h = 0.73. 



LDC, HDC, UZC-SSRS2 and CfAN catalogs, and the 
values only diverge at ~90 Mpc where the linking length 
becomes large enough that the selection biases of the par- 
ent surveys will have a significant effect on the identified 
groups. The number density of groups in the SDSSDR3 
and 2PIGG catalogs is approximately constant over the 
entire range of distances considered due to the lower flux 
limits of the surveys. We note that the HDC catalog 
contains fewer groups than the LDC catalog beyond ^30 
Mpc, as expected. The observed peak at ~70 Mpc cor- 
responds to the location of the great attractor, and is 
enhanced by the distance corrections discussed in i l3.ll 



6. SUMMARY 

We have presented two catalogs of groups in the 
2MASS Redshift Survey, identified using a variable- 
linking- length percolation algorithm (HG82). We dis- 
cussed the effect of the variation of the input parameters, 
Do and Vq, on t he number of groups obtained. As demon- 
strated in H5.4I above, we see that the correct choice of 
parameters depends on the purpose of the catalog, and 
full phase-space information of each galaxy is required 
to understand the most suitable choice of parameters to 
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find virialized groups. 

We justify the choice of two pairs of parameters: (Do, 
Vo) = (1.63 Mpc, 399 km s -1 ), corresponding to a den- 
sity contrast 5p/p = 12 (LDC catalog), and (Dq, Vq) = 
(0.89 Mpc, 350 km s _1 ), corresponding to a density con- 
strast Sp/p = 80 (HDC catalog). We show that the latter 
choice of parameters identifies the largest nearby clusters 
individually, while many of these groups are merged with 
the former parameter choice. 

We compute virial and projected mass estimates for 
the clusters under the assumptions that the identified 
groups have spherical symmetry and that the light traces 
the distribution of the mass. We find that the projected 
mass estimates give mass functions in agreement at the 
3cr-level with Press-Schechter theory, although given the 
dependence of the algorithm on Vo, lcr-agreement is not 
expected. 

We calculate corrected K-band luminosities for each 
cluster, and use these to estimate the mass-to-light ratios 
and corresponding values of Q m . The values predicted 
using the virial mass estimator in both the LDC and 
HDC catalogs are significantly sm aller than the 3-year 
WMAP result l|Snergel et al.ll20f)l of Q m = 0.238^;^ 
motivating the use of the projected mass estimator over 
the virial mass estimator in subsequent analysis. The 
projected mass estimates of groups in the LDC catalog 
produce to a value of fl m = 0.229^o'qi2i which agrees 
with the WMAP result at the la-level. The HDC cat- 
alog significantly under-predicts the WMAP value, sug- 
gesting that by only including groups with density con- 
trasts Sp/p > 80 we are underestimating the total mass 
in groups. 

The distribution of velocity dispersions of groups in 
the 2MRS LDC catalog is i n agreement with the groups 
in the UZC-SSRS2 catalog <|Ramella et al.ll2002j) as well 
as CfAN group catalog. The 2MRS HDC catalog veloc- 
ity dispersions are in agreement with the groups of the 
CfAN catalog, however we find that there is a statisti- 
cally significant difference between the distributions of 
velocity dispersions in the 2MRS HDC and UZC-SSRS2 
catalogs. 

The group catalogs presented in this paper provide an 

17 See follow-up work, Crook et al. (2007a), In preparation. 

18 See Crook et al. (2007b), In preparation. 



estimate of cluster locations and masses without the ne- 
cessity to assume an intrinsic mass-to-light ratio. We 
use and discuss the results of the clustering analysis in 
follow-up work, 18 including the development of a flow- 
field model to estimate the discrcpcncy between the ex- 
pected flow of the local group and t he observed dipole in 
the c osmic microwave background (Be nnett et al.lfl~996l 
120031) . These group catalogs form the basis for a map of 
baryonic density enhancements in the nearby Universe, 
which can be compared with flow-field maps developed 
using peculiar-velocity surveys in order to infer the pres- 
ence and location of dark matter that is not correlated 
with luminous matter. 
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GROUP CATALOGS 



TABLE Al 

Low-density-contrast (LDC) catalog of groups in the 2MASS Redshift Survey 
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NOTE. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. This catalog has been produced using parameters (Dq, Vq) 
= (1.63 Mpc, 399 km s~ 1 ), corresponding to the density contrast 8p/p = 12. We assume h = 0.73 where a value is required. 

a No. of group members (including those generated from the population of the plane). The number derived from the galactic-plane population is contained in parentheses. 

Mean (corrected) group distance. 
C Mcan heliocentric group velocity. 
"^Line-of-sight velocity dispersion. 
e Projected virial radius, 
f 

Log of the virial mass in solar units. 
^Log of the projected mass in solar units. 

1 Log of the (projected) mass-to-light ratio in solar units. 

J Semi-major axis of the ellipse fit to the group at the 75*^— percentile level (measured in arcminutes). 

Axis-ratio of ellipse fit to the group members. 
^Position angle of semi-major axis of ellipse fit to the group members; measured from north toward east. 
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High-density-contrast (HDC) catalog of groups in the 2MASS Redshift Survey 
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Mean (corrected) group distance. 
Mean heliocentric group velocity. 
Line-of-sight velocity dispersion. 
Projected virial radius. 
Log of the virial mass in solar units. 
^Log of the projected mass in solar units. 

1 Log of the (projected) mass-to-light ratio in solar units. 

^ Semi-major axis of the ellipse fit to the group at the 75*^— percentile level (measured in arcminutes). 
Axis-ratio of ellipse fit to the group members. 

Position angle of semi-major axis of ellipse fit to the group members; measured from north toward east. 
Q Group number from LDC catalog that encompasses all members of this group. 
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TABLE A3 

Groups in the LDC catalog with 50 or more members 
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76 (0) 


84.62 


6143 


532.0 


4.23 


15.118 


15.445 


2.255 




127 


01 h 57 m 18 5 s 


+34°19'08" 


71 (0) 


66.22 


4905 


388.1 


3.05 


14.702 


14.966 


2.003 


Perseus-Pisces (A262) 


c 
d 


811 


ll' , 45 m 24.1 s 


+20° 08' 13" 


56 (0) 


95.74 


6548 


550.8 


1.83 


14.785 


14.825 


1.691 


Coma (A 1367) 


Xi 

aa 


292 


04 h 14 m 23.7 s 


+36°58'31" 


51 (0) 


81.93 


6017 


286.1 


3.49 


14.496 


14.810 


1.826 




O 
nr, 



NOTE. — The LDC group catalog was created using (Dq, Vq ) = (1.63 Mpc, 399 km s corresponding to a density contrast Sp/p = 12. Wc assume h = 0.73 where a value is required. 

a No. of group members (including those generated from the population of the plane). The number derived from the galactic-plane population is contained in parentheses. 



Mean (corrected) group distance. 
Mean heliocentric group velocity. 

f-sight velocity dispersion. 
Projected virial radius. 
Log of the virial mass in solar units. 
Log of the projected mass in solar units. 
Log of the (projected) mass-to-light rati- 
Composition of group based on known g, 
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TABLE A4 

Groups in the HDC catalog with 25 or more members 
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a 
it 


"R A 


Dec 


a 

Members 


b 

Distance 


VG 


<Jp 


e 

Hp 


1™ r m v 1 f 


, r m p 1 g 


[ M P /L ]i 

log [^7^J 


Corresponding 


Identified. 












(Mpc) 


(km/s) 


(km/s) 


(Mpc) 


Group #J 
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O 


141 


02 /l 25 m 24.5 s 


+42°05'38" 


29 (0) 


76.52 


5642 


563.1 


1.36 


14.675 


14.801 


2.170 


229 


Perseus- Pisces (A347) 


78 


01^26 m 33.5 s 


-01°34'42" 


29 (0) 


70.73 


5293 


442.7 


0.98 


14.321 


14.500 


1.940 


95 


A194 


O 


56 


01 /i 09 m 52.6 s 


+32°43'06" 


29 (0) 


68.38 


5064 


368.8 


1.30 


14.286 


14.442 


1.808 


79 




c 



807 


13^48 m 47.9 s 


-30°18 / 54 // 


28 (0) 


64.59 


4596 


432.8 


1.20 


14.391 


14.387 


1.751 


956 


Centaurus (A3574) 




763 


13^17 m 41.1 s 


-16°36'08" 


28 (0) 


33.44 


2332 


496.4 


1.09 


14.469 


14.900 


2.759 


852 






432 


07^21 m 47.7 s 


-30°04'19" 


28 (4) 


28.69 


2054 


204.9 


2.19 


14.004 


14.221 


2.247 


391 


NGC2280 Cluster 




282 


04 /i 22 m 27.8 s 


+36°43'26" 


28 (0) 


82.20 


6032 


266.6 


1.79 


14.145 


14.255 


1.520 


292 






732 


12 /i 56 m 09.5 s 


-13°37'44" 


27 (0) 


65.56 


4510 


319.2 


1.34 


14.175 


14.529 


1.926 


883 






295 


04 /i 35 m 49.4 s 


-58°57'46" 


26 (0) 


15.84 


1226 


188.9 


1.50 


13.769 


14.019 


2.077 


391 


Dorado 





Note. — The HDC group catalog was created using (Dq, Vo ) — (0.89 Mpc, 350 km s ), corresponding to a density contrast Sp/p = 80. We assume h = 0.73 where a value is required. 

a No. of group members (including those generated from the population of the plane). The number derived from the galactic-plane population is contained in parentheses. 
^Mean (corrected) group distance. 

Mean heliocentric group velocity. 
^Linc-of-sight velocity dispersion. 

Projected virial radius. 

f 

Log of the virial mass in solar units. 
^Log of the projected mass in solar units. 

1 Log of the (projected) mass-to-light ratio in solar units. 

J Group number from LDC catalog that encompasses all members of this group. 
Composition of group based on known galaxy clusters and supcrclusters . 
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TABLE A 5 

Groups in the LDC catalog and their members 



Name 


RA 


Dec 




m K b 


Distance 


Corresponding 






(km/s) 




(Mpc) 


Group # d 


Group 1 


000009.14+3244182 


00"00 m 09.0 s 


+32°44'18" 


10372 


10.61 


140.48 


1 


000028.80+3246563 


00 ,I 00 m 28.8 s 


+32°46'56" 


9803 


11.09 


132.77 


1 


235950.52+3242086 


23 h 59 m 50.5 s 


+32°42'09" 


10086 


11.12 


136.61 


1 


Group 2 


000104.78+0432261 


00 h 01 m 04.7 s 


+04° 32' 26" 


9151 


11.09 


123.42 


None 


000221.65+0405230 


00 h 02 m 21.6 s 


+04°05'23" 


8606 


11.10 


116.03 


None 


235858.87+0338045 


23 h 58 m 58.8 s 


+03°38'04" 


8854 


10.93 


119.43 


None 


Group 3 


000037.94+2823041 


00 fe 00 m 37.8 s 


+28°23'04" 


8705 


10.46 


117.75 


1258 


000046.96+2824071 


00 h 00 m 46 gs 


+28°24'07" 


8764 


10.41 


118.55 


1258 


000433.73+2818059 


Oo' l 04 m 33.6 s 


+28°18'06" 


8785 


10.62 


118.79 


None 


235828.41+2802025 


23' , 58 m 28.3 s 


+28°02'03" 


9145 


10.94 


123.72 


1258 


235943.72+2817251 


23 h 59 m 43.6 s 


+28°17'25" 


9073 


10.71 


122.74 


1258 


Group 4 


000100.43+0614312 


00 fe 01 m 00.3 s 


+06°14'31" 


5324 


10.03 


71.77 


None 


000348.85+0728429 


00 ,l 03 m 48.9 s 


+07°28'43" 


5241 


9.67 


70.63 


None 


000649.47+0837425 


00 ,I 06 m 49.5 s 


+08°37'42" 


5257 


10.58 


70.82 


None 


235651.54+0530303 


23 h 56 m 51.6 s 


+05°30'30" 


5405 


11.05 


72.90 


None 


Group 5 


000105.97-5359303 


00"01 m 06.0 s 


-53°59'30" 


9423 


10.57 


129.47 


None 


000310.64-5444562 


00 ,I 03 m 10.6 s 


-54°44'56" 


9767 


10.34 


134.20 


None 


000311.27-5444588 


00 h 03 m 11.3 s 


_54°44'59" 


9767 


10.35 


134.20 


None 


Group 6 


000329.22+2721063 


00 fe 03 m 29.1 s 


+27°21'06" 


7690 


11.02 


103.97 


4 


000548.43+2726579 


00 ,i 05 m 48.3 s 


+27°26'58" 


7531 


10.95 


101.81 


4 


000717.10+2740421 


00 ,I 07 m 17.1 s 


+27°40'42" 


7550 


11.13 


102.05 


4 



Note. — The complete version of this table is in the electronic edition of the Journal. The table in the printed edition is only intended 
as a guide to its content. The LDC group catalog was created using (Do, Vq) = (1.63 Mpc, 399 km s~ x ), corresponding to a density 
contrast Sp/p = 12. 

a Heliocentric velocity. 

Corrected distance, assuming h = 0.73. 
C Apparcnt K magnitude. 

Corresponding group number assigned to this galaxy when in the HDC catalog. 
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TABLE A6 

Groups in the HDC catalog and their members 



Name 


RA 


Dec 




m K b 


Distance 


Corresponding 






(km/s) 




(Mpc) 


Group # d 


Group 1 


000009.14+3244182 


00"00 m 09.0 s 


+32°44'18" 


10372 


10.61 


140.48 


1 


000028.80+3246563 


00 ,I 00 m 28.8 s 


+32°46'56" 


9803 


11.09 


132.77 


1 


235950.52+3242086 


23 h 59 m 50.5 s 


+32°42'09" 


10086 


11.12 


136.61 


1 


Group 2 


000001.68+4716282 


00"00 m 01.7 s 


+47° 16' 28" 


5017 


9.68 


68.91 


11 


000012.95+4657543 


00 h QO m 13 qs 


+46°57'54" 


5366 


10.84 


73.57 


11 


000426.65+4729250 


Oo' i 04 m 26 .6" 


+47° 29' 25" 


5269 


10.52 


72.27 


11 


000527.96+4632371 


00 ,i 05 m 28.0 s 


+46°32'37" 


4971 


11.02 


68.21 


11 


000723.79+4702265 


00' , 07 m 23.8 s 


+47° 02' 27" 


5313 


9.93 


72.81 


11 


000724.58+4659195 


00 h 07 m 24.6 s 


+46°59'20" 


5097 


10.87 


69.91 


11 


235247.40+4648138 


23 ,l 52 m 47.3 s 


+46°48'14" 


5047 


10.65 


69.35 


11 


235401.14+4729225 


23 h 54 m 01 1S 


+47°29'22" 


5202 


10.39 


71.45 


11 


235526.14+4716485 


23 h 55 m 26.1 s 


+47°16'49" 


5348 


10.98 


73.38 


11 


235915.79+4653213 


23 h 59 m 15.8 s 


+46°53'21" 


5021 


9.46 


68.95 


11 


Group 3 


000126.77+3126016 


00"01 m 26.7 s 


+31°26'02" 


4948 


10.23 


67.25 


9 


000130.05+3126306 


00 ,l 01 m 30.0 s 


+31°26'31" 


4767 


10.47 


64.83 


9 


000308.87+3102108 


00 ,i 03 m 08.9 s 


+31°02'11" 


4797 


11.16 


65.20 


9 


000424.49+3128193 


00 h 04 m 24.5 s 


+31°28'19" 


4958 


11.14 


67.36 


9 


Group 4 


000329.22+2721063 


00 h 03 m 29.1 s 


+27°21'06" 


7690 


11.02 


103.97 


6 


000548.43+2726579 


00 ,i 05 m 48.3 s 


+27°26'58" 


7531 


10.95 


101.81 


6 


000717.10+2740421 


00 h 07 m 17.1 s 


+27°40'42" 


7550 


11.13 


102.05 


6 


Group 5 


000457.78+0507245 


00"04 m 57.8 s 


+05°07'24" 


5357 


11.19 


72.17 


7 


000527.66+0513204 


00 ,l 05 m 27.6 s 


+05°13'20" 


5294 


10.08 


71.31 


7 


000640.35+0506483 


00 ,1 06 m 40.3 s 


+05°06'48" 


5371 


11.22 


72.34 


7 



Note. — The complete version of this table is in the electronic edition of the Journal. The table in the printed edition is only intended 
as a guide to its content. The HDC group catalog was created using (Do, Vq) = (0.89 Mpc, 350 km s _1 ), corresponding to a density 
contrast Sp/p = 80. 

a Hclioccntric velocity. 

^Corrected distance, assuming h = 0.73. 
C Apparent K magnitude. 

^Corresponding group number assigned to this galaxy in the LDC catalog. 
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